Pis'ma v ZhETF 

Quark cluster contribution to cumulative proton emission in 
fragmentation of carbon ions. 

B. M. Abramov, P. N. Alekseev, Yu. A. Borodin, S. A. Bulychjov, I. A. Dukhovskoy, A. I. Khanov, A. P. Krutenkova 1 ^, 
V. V. Kulikov, M. A. Martemyanov, M. A. Matsyuk, E. N. Turdakina 

Institute for Theoretical and Experimental Physics, 117218 Moscow, Russia 



CO 

o 

CN 

< 

CO 
CN 



X 
=5 



> 

o 

CN 
CN 

O 

CO 



X 



In the FRAGM experiment at heavy ion accelerator complex TWAC-ITEP, the proton yields at an angle 
3.5° have been measured at fragmentation of carbon ions at To = 0.6, 0.95 and 2.0 GeV/nucleon on beryllium 
target. The data are presented as invariant proton yields on cumulative variable x in the range 0.9 < x < 
2.4. Proton spectra cover six orders of invariant cross section magnitude. They have been analyzed in the 
framework of quark cluster fragmentation model. Fragmentation functions of quark-gluon string model are 
used. The probabilities of the existence of multi-quark clusters in carbon nuclei are estimated to be 8-12% for 
six-quark clusters and 0.2-0.6% for nine-quark clusters. 



Introduction. Since the discovery of cumulative 
effect in seventies [TJ [2J, the question of the nature 
of cumulative particles is still under discussion. 
Cumulative particles are produced in interactions with 
nuclei in kinematic region forbidden for interaction 
with free nucleon. Few hypothesis of their origin have 
been considered. Among them there are fluctuations of 
nuclear matter [3], clusters [I], few- nucleon correlations 
[5], excited dibaryons [5], etc. It was pointed out [71 [5] 
|that multiple scattering in nuclear matter can also 
contribute. Modern theoretical approaches (see, e.g., 
[9]) connect the cumulative effect with contribution 
from multi-quark states. These states (first of all, 
six-quark states) are considered for nuclear matter 
phase transitions at high densities. At intermediate 
energies |10| . the existence of multi-quark clusters 
in cold and hot baryonic matter blurs the boundary 
between hadronic and quark-gluon phases. Probability 
of the two-nucleon fluctuation in 12 C nucleus was 
for the first time estimated in [TT] based on the 
experimental data on cumulative pion production. It 
is in qualitative agreement with theoretical prediction, 
which was obtained later within quark cluster model 
|12| . Similar approach, in which fragmentation functions 
of quark clusters were calculated within the quark- 
gluon string model, was successfully used in [3] for the 
description of K~ , n- and antiproton inclusive spectra 
in hadron-nucleus interactions. 

Ion beams open new ways for the cumulative effect 
study. One of the authors of [3], A.B. Kaidalov, 
suggested to use the data from our experiment FRAGM 
[inilll] for the experimental estimation of the admixture 
of the multi-quark cluster state in nuclear matter. In this 
experiment, cumulative protons are measured in inverse 
kinematics, i.e. in the fragmentation region of projectile 
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nucleus. Such measurement has definite advantages over 
measurements in the target fragmentation region. First, 
relativistic boost of forward going protons increases 
considerably proton detection solid angle in the rest- 
frame of fragmented nucleus at fixed acceptance of 
a detection system in laboratory frame. Second, in 
the inverse kinematics there are no problems with a 
detection of protons which are at rest in projectile 
rest frame because they have lab momentum close 
to momentum per nucleon of the projectile. This 
allows to clearly observe nucleon-nucleon component 
of the projectile and use it for normalization that is 
impossible for target fragmentation. However inverse 
kinematic creates additional problems, connected with 
measurements of the protons with momentum few times 
higher than momentum per nucleon of the projectile. 

Experiment. In the FRAGM experiment at the 
accelerator complex TWAC (Tera-Watt Accumulator) 
ITEP, the yields of nuclear fragments are studied in 
carbon fragmentation on beryllium target: 



2 C + Bc^f+X. 



(1) 



The main goal of the experiment is to collect data at 
high momenta of nuclear fragments f. In this paper the 
proton spectra from the reaction (1), obtained at carbon 
kinetic energies To = 0.6, 0.95 and 2.0 GeV/nucleon 
are analyzed. A measurement of momentum spectra 
at different energies gives a possibility to study 
energy dependence of their parameters. Experimental 
set-up (Fig. 1) is comprised of two-step magneto- 
optical channel, placed at an angle 9 = 3.5° with 
respect to internal ion beam of the ITEP accelerator. 
Narrow vertical strip of 200 u.m Be-foil was used as a 
target which allows to have simultaneously both high 
luminosity due to multiple passage of ions through the 
target and small sizes of the source for full usage of 
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beryllium, boron and carbon isotopes are clearly 
seen. At each projectile energy the proton yields are 
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Fig.l. Set-up of the FRAGM experiment. TWAC MR 
is the TWAC main ring, M is a monitor, Q1-Q5 are 
quadrupoles, BM1 and BM2 are bending magnets, Col 
and Co2 are collimators, CF1, CF2, HI, C2 and C3 are 
scintillator counters. 



high momentum resolution of the channel. The first 
step of the channel consists of quadrupole doublet Ql, 
Q2, bending magnet BM1 and achromatic correction 
quadrupole Q3, placed in the first focus of the channel. 
The second step contains bending magnet BM2 and 
quadrupole doublet Q4, Q5. This step refocuses the 
target image from the first focus to the second one, 
sixteen meters downstream. Scintillator hodoscope of 
twenty vertical and eight horizontal elements with 
sizes 20x1x1 cm 3 for beam profile measurements is 
placed in the first focus. This hodoscope allows to 
improve fragnrent momentum resolution up to 0.2%, 
using focusing properties of magneto-optical channel 
with full momentum acceptance of ±3%. Scintillator 
counters CF1, CF2 and C2, C3 for amplitude and 
time-of-flight (TOF) measurements are placed in the 
first and in the second focuses. Each counter is seen by 
two PMT's from the opposite sides to use their signals 
in mean-time determination. Coincidence of signals 
from the counters of the first and the second focuses 
gives a trigger, which starts the data transfer from 
CAM AC system to a computer under LINUX. Readout 
software is based on ROOT [T5] package. As a monitor 
M, a telescope of three scintillator counters, which view 
the target at an angle of about 2° is used. 

Protons are selected on correlation plot of signal 
charge from a counter (a function of fragnrent charge) 
vs TOF (a function of fragment mass). As an example, 
correlation plot for projectile energy of 0.6 GeV and 
momentum of magneto-optical channel of 2.15 GeV/c 
is shown in Fig. 2. Groups of events corresponding 
to outgoing fragments: hydrogen, helium, lithium, 
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Fig. 2 Correlation plot of signal charge from one of 
scintillation counters (QDC) vs time-of-flight (TDC) at 
a base of 16 m. Projectile energy is 600 MeV/nucleon, 
fragment rigidity is 2.15 GeV/c. 

measured by momentum scan of the channel with a step 
25-200 MeV c _1 . Set-up efficiency is calculated with 
GEANT4-based [T^ simulation program. The program 
traces particles in the magneto-optical channel taking 
into account multiple scattering effects, ionization 
losses and absorption in the detector materials. 

Data analysis. The invariant cross sections of 
proton yield <Ji nv = (E / po)d 2 a / dxd(p 2 ) as a function of 
cumulative variable x — p/po for projectile energies To 
=0.6, 0.95 and 2.0 GeV/nucleon are shown in Fig. 3 - 5. 
Here po is projectile momentum per nucleon, p is proton 
momentum in the laboratory frame, pt is its transverse 
component with respect to projectile. The data cover six 
orders of invariant cross section magnitude. It is three 
orders of magnitude more than in the most sensitive 
previous experiment |17) at 1.05 GeV/nucleon. In the 
region of the maximum (x ~ 1.0) the shapes of the 
spectra are close to Gaussians as predicted by statistical 
models. However, already at x > 1.3 the spectra become 
exponentials, which is typical for cumulative processes. 

As was already mentioned, the most successful 
approach to the problem of cunrulative particle 
production is the quark cluster model [2]. This model 
was used to describe yields of cunrulative pions, kaons 
and antiprotons. However, for protons such analysis was 
not performed. Reliable identification of one-nucleon 
component achievable in inverse kinematics allows to 
apply this analysis also for protons. In the framework 
of this model clusters existing in a nucleus, consist 
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Fig. 3. Points are invariant cross section (Ji nv of proton 
yield at projectile energy of 0.6 GeV/nucleons in 
arbitrary units as function of x variable. Solid curve is 
a fit to the data in quark-cluster model. Contributions 
from one-, two- and three-nucleon clusters are shown 
with dashed, dotted and dash-dotted lines, respectively. 
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Fig.4. The same as Fig. 3, but at 0.95 GeV/nucleon. 



of 3fc (k = 1, 2, 3,...) valent quarks. Conventional 
nucleon component of the nucleus corresponds to k = 
1. The probabilities Wk of such clusters in a nucleus 

A-\ 

with A nucleons are normalized by Y] Wk — 1. As 

k=l 

cluster contribution into observed processes falls with 
increasing fc, we limit ourselves with k = 1, 2, 3 only 
and represent the invariant cross section as a sum of 
three components: 



a inv oc Gwig(x,p 2 ) + w 2 b 2 (x,p 2 ) + w 3 b 3 (x,p 2 ), (2) 
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Fig. 5. The same as Fig. 3, but at 2.0 GeV/nucleon. 

where functions <?, b 2 and 63 are fragmentation functions 
of quark clusters into protons, g is a gaussian 



g{x,p 2 t ) = 

ex P (-0.5((l - A) - xfjol) exp(-0.5p 2 /a 2 ), 



(3) 



&2 and &3 are calculated in the framework of the quark- 
gluon string model: 




B 2 (x/2) 3 {l-x/2) 3 exp(-a lP 2 t ), 

0, x^[0,2], 

B 3 (x/3f(l-x/3feM~a 2 p 2 ), 



(4) 



(5) 



0, x(f. [0,3]. 

The values G, B 2 and B 3 are known normalization 
constants; the first one is defined by Gaussian 
normalization 

G = (4V27rcr a; CTp) _1 , cr p = cr x m p p /(T + m p ), (6) 

and two others, B 2 and B 3 , can be found from: 



bi(x,p 2 )dxdp 2 = z/2, i = 2,3. 



(7) 







Transverse momentum dependence is not predicted 
within quark gluon string model. We have got the 
invariant cross section slopes on p 2 from the data [17] 
and extrapolated them to the region of our experiment, 
the values ai ~ 5 GeV _2 c 2 and a 2 = 3 GeV~ 2 c 2 were 
used. 

The results of the data fit with the formula (2) at 
projectile energies 0.6, 0.95 and 2.0 GeV are shown 
in Fig. 3-5. The fitted curve is shown by solid line. 
The contributions of one-nucleon (3q) component, two- 
nucleon (6q) and three-nucleon (9q) clusters are given 



by dashed, dotted and dashed-dotted lines, respectively. 
The fit parameters are mean value (1-A) and r.m.s. 
of Gaussian (a x ), and also probabilities w 2 and W3, 
connected to w\ with the relation w\ + u>2 + W3 
= 1. The obtained probabilities W2 and W3 for carbon 
nucleus are given in the Table. The two-nucleon cluster 

Table. Results from the fits of Fig. 3 - 5. To and po 
are kinetic energy and momentum of the projectile per 
nucleoli, x m ax is maximal value of x, reached in this 
experiment, probabilities u>2 and W3 for carbon nucleus 
are defined in the text. Statistical errors of the fit are 
given in parentheses. 



To, 
GeV/n 


Vo, 

GeV/c/n 


•Emax 


w 2 


w 3 


0.6 


1.22 


1.95 


.077(10) 


.004(2) 


.95 


1.6 


2.4 


.119(17) 


.002(1) 


2.0 


2.72 


2.15 


.098(18) 


.006(1) 



probability, estimated at different projectile energies 
varies within 7.7 - 11.9%, while the three-nucleon one 
is within 0.2 - 0.6%. They are compatible both with 
given statistical errors of the fit and with expected 
independence of these probabilities on projectile energy. 
It should be emphasized that obtained probabilities for 
carbon nucleus could be considered only as estimates 
because of difficulties to take into account systematic 
uncertainties of the theoretical approach. First of all 
the fragmentation functions (4) and (5) were justified 
in the quark-gluon string model only in the boundary 
regions near x=0, x=2 and x=3. The calculation of 
fragmentation functions in the whole x range is still 
unsolved problem of the model. The fragmentation 
functions used here are the simplest (and widely used) 
ones, satisfying above mentioned boundary conditions. 
Moreover, in Fig. 3 - 5 at x > 1.3 one can see wave- 
like behavior of the resulting fit curve while the data 
demonstrate purely exponential fall off with x. 

This can be considered as indication that the model 
does not take into account some small effects which 
smear fragmentation functions behavior. One of them 
is evident, it is internal motion of two- and three- 
nucleon clusters in a nucleus which has not been taken 
into account in this approach. It can be expected that 
improvement of multi-quark cluster approach will allow 
to overcome these difficulties. The important result of 
the presented analysis is the demonstration of possibility 
to use data on cumulative proton yield in the inverse 
kinematics for estimates of these probabilities. Obtained 
value of W2 is close to the value of 6%, obtained in 
|11) from cumulative pion production and to theoretical 



predictions of 12.5% given in [12] . However, the value 
of W3 is much smaller than 2.6% predicted in |12| . The 
values u>2 (103) are not far from probabilities of two- 
nucleon (three-nucleon) correlations in nuclei, obtained 
in the TJNAF experiment on A(e,e') inclusive electron 
scattering on nuclei [H], (19.3 ± 4.1)% and (0.55 ± 
0.17)% for carbon nucleus. Reasonable agreement of 
the results of these experiments can be considered as 
evidence for unique nature of quark clusters and short- 
range nucleon correlations in nuclei. 
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